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ABSTRACT 

We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet 
to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the 
Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple 
wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of 
all galaxies for which FIR Herschel SPIRE observations are publicly available. Our SEDs combine the 
Herschel photometry with mult i- wavelength data from Spitzer^ GALEX^ Swift UVOT, and 2MASS. 
While the shapes of the SEDs are broadly similar across our sample, strongly interacting galaxies 
typically have more mid-infrared emission relative to their near-infrared and FIR emission than weakly 
or moderately interacting galaxies. We modeled the full SEDs to derive host galaxy star formation 
rates (SFR), specific star formation rates (sSFR), stellar masses, dust temperatures, dust luminosities, 
and dust masses. We find increases in the dust luminosity and mass, SFR, and cold (15-25 K) dust 
temperature as the interaction progresses from moderately to strongly interacting and between non- 
interacting and strongly interacting galaxies. We also find increases in the SFR between weakly and 
strongly interacting galaxies. In contrast, the sSFR remains unchanged across all the interaction 
stages. The ultraviolet photometry is crucial for constraining the age of the stellar population and 
the SFR, while dust mass is primarily determined by SPIRE photometry. The SFR derived from the 
SED modeling agrees well with rates estimated by proportionality relations that depend on infrared 
emission. 

Subject headings: infrared: galaxies, galaxies: interactions, galaxies: photometry, galaxies: star for- 
mation, ultraviolet: galaxies 



1. INTRODUCTION 

Galaxy evolution is believed to be heavily infiuenced 
by interactions between galaxies, both for local systems 
and for distant objects at earlier cosmological times. In 
the canonical view, interactions between galaxies have 
three primary observable effects. In the most dramatic 
cases, interactions stimulate star formation in a burst of 
activity that is presumed to power the high infrared (IR) 
luminosities typically seen in such systems. Many local 
ultra-luminous IR galaxies (L > 10^^ Lq; ULIRGs) and 
luminous IR galaxies (10^^ L© < L < 10^^ L©; LIRGs) 
show evidence of galaxy interactions (e.g., Veilleux et 
al. 2002). Similarly, their high-redshift counterparts, 
sub-millimeter galaxies, first detected by SCUBA and 
now studied extensively by the Spectral and Photometric 
Imaging Receiver (SPIRE) on the Herschel Space Obser- 
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vatory^^ are thought to be predominantly mergers (e.g., 
Blain et al. 1999), although the relative contribution of 
mergers of different stages to their numbers is still an 
open question (e.g., Hayward et al. 2012a, 2012b). 

The second effect is that interactions significantly af- 
fect the subsequent evolution of galaxies, which may lead 
to significant changes in their morphology. Disturbed 
galaxies have long been associated with mergers (e.g., 
Toomre & Toomre 1972). Numerical simulations of in- 
teractions (e.g., Hopkins et al. 2006; Hopkins 2012; Mi- 
hos & Hernquist 1994, 1996; Barnes 1992; Barnes & 
Hernquist 1996; Sanders 1999) show a variety of mor- 
phological distortions as well as variable amounts of star 
formation. The simulations also demonstrate the com- 
plexity of the problem: the degree of induced activity 
and distortion varies greatly with the parameters of the 
encounter, the phase of the interactions, the molecular 
gas content ( "wetness" ) , and the mass of the progenitor 
galaxies among many other properties. 

Third, the canonical picture, as seen in many simula- 
tions (e.g., Di Matteo et al. 2005, Springel et al. 2005), 
involves merger-driven gas infiow to the central regions, 
resulting in heightened activity of the central supermas- 
sive black hole as well as starburst activity due to the in- 
creased central gas density and possibly turbulence. The 
process in principle converts a low- luminosity nucleus 
into an active galactic nucleus (AGN) but one whose lu- 

^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA. 



2 



Lanz et al. 



minosity might range widely depending on the stage of 
the interaction. Indeed, observations of merging galaxies 
over the years have tended to provide evidence support- 
ing the conclusion that, at least on a statistical level, 
interactions trigger an enhancement in the formation of 
stars as well as nuclear activity. However, the recent lit- 
erature includes works that argue both for and against a 
strong connection between nuclear activity and mergers 
(e.g., Li et al 2008, Kocevski et al 2012, Ellison et al 
2011, Scudder et al 2012, Silverman et al 2011). 

Therefore, despite many previous studies (e.g., Dahari 
1985, Sanders et al 1988, Kewley et al 2001, Lambas et 
al 2003), both observational and through simulations, 
our understanding of the evolution of the physical ac- 
tivity during the course of a galaxy-galaxy interaction 
remains incomplete. In the past decade, two new devel- 
opments have dramatically changed our understanding 
of star formation and accretion activity around galactic 
nuclei, which are the two dominant processes at work in 
controlling the observed emission. The first is the suc- 
cess of space missions, in particular, the Spitzer Space 
Telescope^ (Werner et al 2004) and the Herschel Space 
Observatory (Pilbratt et al 2010) in the IR, as well as 
the Galaxy Evolution Explorer^^ [GALEX] Martin et 
al 2005) and Swift (Gehrels et al 2004) in the ultraviolet 
(UV) , providing photometry across the complete spectral 
range from UV to millimeter wavelengths. Most impor- 
tantly, detailed imaging and high sensitivity photometry 
now available at the critical far-infrared (FIR) emission 
peak resulting from warm, luminous dust heated by star- 
bursts provides crucial information regarding dust heat- 
ing and embedded star formation. The combination of 
UV and IR observations is essential to obtain a complete 
census of recent and ongoing star formation by capturing 
both the unobscured and obscured emission from young 
stars. 

The second development has been the success of com- 
putational codes. We have new tools for the derivation of 
galaxy properties including masses, star formation rates, 
and interstellar medium (ISM) parameters from global 
fits, which allow self-consistent measurements of criti- 
cal parameters combining stellar evolution models (e.g., 
Bruzual & Chariot 2003) with radiative transfer through 
a dusty ISM (e.g.. Chariot & Fah 2000). A second set 
of tools uses sophisticated hydrodynamic computational 
codes to simulate interactions (e.g., GADGET - Springel 
2005), while simultaneously new radiative transfer mod- 
els can compute the predicted emission from these evolv- 
ing interacting systems (e.g., SUNRISE - Jonsson 2006). 

It is important to recognize that observational biases 
can be significant. Due to the long timescales of an in- 
teraction (typically 10^~^ years), observers rely on stud- 
ies of a range of interacting systems to reconstruct a 
likely sequence of events. Moreover, determining the ex- 
act phase of any particular observed interaction from its 
morphology is uncertain at best, because the appearance 
of a system at a given interaction phase also depends 
on the specific geometry of the encounter, the masses 
of the galaxies, metallicity, molecular gas content, and 

^ Spitzer is operated by the Jet Propulsion Laboratory, Califor- 
nia Institute of Technology under a contract with NASA. 

GALEX is operated for NASA by the California Institute of 
Technology under NASA contract NAS5-98034. 



(not least) previous interaction histories (e.g., Di Mat- 
teo et al 2007). Systems are ordered into an evolution- 
ary sequence using intuition provided by simulations and 
physical models, which are themselves based on observa- 
tions of particular systems. Selection criteria, however, 
can introduce a bias for more luminous, morphologically 
disturbed systems and, hence, towards the most active 
phases of interactions. Therefore, a selection criterion 
not associated with either morphological disturbance or 
degree of activity is critical for obtaining a sample con- 
taining systems throughout the full interaction sequence. 

We have undertaken a program to take advantage of 
all these developments: full, mult i- wavelength datasets 
of an interacting galaxy sample selected with few biases; 
hydrodynamic simulations; and radiative transfer model- 
ing, in a systematic effort to better understand systems 
across a range of interaction stages and to iterate re- 
finements to the various modeling and radiative transfer 
codes. We have chosen a representative sample of ob- 
jects spanning the interaction sequence, obtained their 
full spectral energy distributions (SEDs), and are com- 
paring the results against a variety of models - based on 
both templates/stellar evolution/radiative transfer and 
on diagnostic features. 

This first paper of the project presents results and con- 
clusions for a sample of thirty-one interacting galaxies 
in fourteen systems for which there are currently com- 
plete multi-wavelength data that can be used to study 
the variations in their star formation and dust heating. 
This paper is organized as follows. We describe the full 
Spitzer Interacting Galaxy Survey (SIGS) sample in Sec- 
tion 2 and the classification of each of the sources in the 
interaction sequence. Section 3 describes the wide range 
of observational photometry used to construct the SEDs. 
It also describes the issues associated with obtaining reli- 
able photometry from the diverse datasets. In Section 4, 
we model the SEDs of these objects. Section 5 discusses 
the variations seen across the interaction sequence and 
constraints imposed by photometry from different instru- 
ments and compares star formation rates derived using 
the entire SED to those from relations depending on one 
or two wavelengths. In Section 6, we summarize our re- 
sults. 

2. THE SPITZER INTERACTING GALAXY SURVEY (SIGS) 

SAMPLE 

2.1. Sample Description 

The SIGS sample was designed to span the full range of 
galaxy interaction parameters by using a sample selected 
strictly on the basis of interaction probability rather than 
morphology, activity, luminosity, or other derivative in- 
dicators. The catalog includes interactions of all types, 
not just those that give rise to obvious morphological pe- 
culiarities and/or nuclear /starburst activity, thus mini- 
mizing morphological biases so we can address the rela- 
tionships between interactions and activity. A selection 
criterion not dependent on visible signs of tidal inter- 
actions is important because of the dependence of the 
response of interacting galaxies on the relative inclina- 
tions of disks (e.g., Toomre & Toomre 1972; D'Onghia 
et al 2010) and the uncertain distribution of dark mat- 
ter around the galaxies (e.g., Dubinski et al 1996, 1999). 
The SIGS sample was based on the Keel-Kennicutt vis- 
ibly selected catalog of interacting spiral galaxies (Keel 



SIGS Paper I 



3 



TABLE 1 
Sample Description 
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et al. 1985, hereafter K85), which selected galaxies based 
on the local density of nearby neighbors and consists of 
bright spiral galaxies having neighbors with typical pro- 
jected separations of 4-5 effective radii. A criterion based 
on the relative recessional velocities |Av| < 600 km s~^ 
was imposed to exclude non-associated, projected pairs. 

In order to resolve structures on scales of a few hun- 
dred pc, we limited the original sample to sources closer 
than cz < 4000 km s~^. To investigate the effects of 
tidal interaction, we added a complementary set to the 
prime sample: the K85 "Arp Sample" with the same 
maximum distance as the Keel-Kennicutt complete sam- 
ple. This set is based on the Arp catalog of peculiar 
galaxies from which K85 selected all objects showing ev- 
idence of tidal interaction not strong enough to disrupt 
the galactic disks (i.e., it does not bias against late stage 
mergers). Although K85 excluded some fainter members 
of the interacting groups (their selection criteria required 
a B band magnitude of Bt < 13.0), we include them in 
order to obtain a complete picture of the activity in the 
different interacting systems. 

The total SIGS sample consists of 103 individual in- 
teracting galaxies in 48 systems. The combined galaxies 
span the range of interaction types, luminosities, and 
galaxy types. SIGS is comprised primarily of spiral- 
spiral interactions, with some spiral-elliptical and spiral- 
irregular interactions. Its set of systems contains both 
major and minor mergers, ranging from systems likely 



to be in first approach (e.g. NGC 3424/NGC 3430) 
through close passages (e.g. M51) to final collision 
(NGC 3690/IC 694), and span a luminosity range from 
1.3x10^^ — 5.1 X 10^^ L0. From this complete sample, 
which has a sufficiently large number of objects to al- 
low us to study statistically the activity in interacting 
galaxies across a wide range of encounter parameters, we 
will be able to study the increase of star formation and 
AGN activity in interacting disk galaxies. As discussed 
in section 1, while there have been a significant number of 
studies probing star formation rate (SFR) enhancement 
and nuclear activity, the importance of the different in- 
teraction parameters in triggering these events is not well 
understood. The SIGS sample provides us with the op- 
portunity to observe a large range of galaxies, including 
very early interaction stages. The level and distribution 
of star-formation in such early stage interactions has not 
been systematically studied before, therefore our sample 
will allow us to identify the initial increase in SFR caused 
by the interaction, as well as identify where this enhance- 
ment is located in the galaxies (i.e. in the central region 
of the galaxy, along the disk, or within tidal features). 
Additionally, the size of our sample also provides us with 
the ability to probe these enhancements for all systems 
as a function of different interaction parameters, such as 
galaxy mass, mass ratio and gas content. A detailed de- 
scription of the SIGS sample along with the analysis of 
the Spitzer data and a presentation of the images and the 
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photometric results is given in Brassington et al. (2013, 
in preparation). 

There are currently fourteen interacting systems from 
the SIGS set which have publicly available observations 
by all the facilities: Herschel (SPIRE and partial cov- 
erage with the Photodetector Array Camera and Spec- 
trometer (PACS)), Spitzer, 2 Micron All Sky Survey 
(2MASS), and either GALEX or Swift, enabling us to 
model their emission from far-UV (FUV) to FIR in 28 
filters when ancillary archival measurements are added. 
Not all galaxies have photometric data in all filters; we 
used as many photometric data as available, generally 
15-25. These galaxies comprise the sample we examine 
in this paper and were selected from the SIGS sample 
on the basis of available SPIRE observations. They are 
listed in Table 1 along with key parameters. 

2.2. Estimating the Interaction Phase 

Toomre & Toomre (1972) were the first to systemat- 
ically model and describe the morphological character- 
istics of interacting galaxies. Using simple simulations, 
they showed that tails and bridges could result from tidal 
forces and reconstructed the orbits that could produce 
the tidal features seen in some of the best known inter- 
acting systems including M51, the Mice (NGC 4676), and 
the Antennae (NGC 4038/4039). Their work also high- 
lighted the close connection between observations and 
modeling: our classification of the interaction stages in 
our sample is based on theoretical descriptions of how 
such interactions are expected to proceed. 

As Rich et al. (2012) have shown, projected distance 
alone is an unreliable indicator of interaction stage. We 
therefore used the Dopita et al. (2002) five-stage scheme 
to classify the interaction stage of our galaxies. By con- 
struction, our sample does not include any Stage 1 galax- 
ies (isolated, non-interacting galaxies). Stage 2 galaxies 
are described as weakly interacting systems, which are 
close on the sky, but show minimal morphological dis- 
tortion. These systems could be either before or after 
the first passage. Stage 3 galaxies, which we call mod- 
erately interacting, show stronger signs of morphological 
distortion and often tidal tails. Depending on the geome- 
try of the encounter, these systems could be in the midst 
of the first or a subsequent passage. Stage 4 (strongly in- 
teracting) galaxies show strong signs of disturbance and 
are therefore in more evolved stage of interaction. Our 
sample falls into these three categories. While the SIGS 
sample has a Stage 5 (coalescence/post-merger systems), 
the sample presented in this paper does not. The SIGS 
sample is roughly equally divided between Stages 2-4, 
while the sample presented in this paper has 7, 14, and 
7. 

This classification method is clearly a statistical 
scheme in the sense that, for each individual galaxy, the 
classification stage does not translate directly to an in- 
teraction phase. However, since the scheme is based on 
morphological appearance of galaxies, it provides a di- 
rect picture of the effect of the interaction on the distri- 
butions of the stellar component of the galaxies and their 
star formation activity. The classification was carried out 
independently for each galaxy in the SIGS sample by six 
collaborators on the basis of appearance alone in Digi- 
tized Sky Survey (DSS) images. Stage 2 galaxies show 
little morphological distortion, while Stage 4 galaxies are 



strongly distorted. Stage 3 galaxies show some distortion 
in the form of tidal features, although their disks remain 
undisturbed. Visible DSS images are best suited for this 
purpose, since they trace on-going star formation as well 
as older stellar populations in a single image. In Ap- 
pendix A, we show representative examples of the galax- 
ies in Stages 2-4. Galaxy groups in which classifications 
differed by more than one stage were re-examined; the 
median of the classifications is used for each galaxy. Ta- 
ble 1 lists the interaction stage for all of the galaxies in 
our sample. 

2.3. Comparison Non- Interacting Sample 

As a comparison sample of non-interacting galaxies, we 
used a subset of the "normal" galaxy sample of Smith 
et al (2007). Smith et al (2007) identified 42 galax- 
ies from the Spitzer Infrared Nearby Galaxies Survey 
(SINGS; Kennicutt et al 2003; Dale et al 2005) of which 
26 were spirals, which had not been subject to strong 
distortions. We were more conservative in our definition 
of non-interacting, by removing galaxies associated with 
clusters or radial- velocity groups, and we removed the 
three that were not observed with SPIRE as part of the 
Key Insights on Nearby Galaxy: a Far Infrared Survey 
with Herschel (KINGFISH; Kennicutt et al 2011). Our 
comparison sample is comprised of 15 galaxies: NGC 
925, NGC 1291, NGC 2841, NGC 3049, NGC3184, NGC 
3521, NGC 3621, NGC 3938, NGC 4236, NGC 4559, 
NGC 4594, NGC 4736, NGC 4826, NGC 5055, and 
NGC 6946. We used the distances provided in Smith 
et al. (2007) and the UV-MIR photometry given in Dale 
et al 2007) and the FIR photometry given in Dale et 
al (2012). 

3. OBSERVATIONS AND DATA REDUCTION 

The sample presented here has a complete set of near- 
infrared (NIR) to FIR photometry observed by 2MASS, 
Spitzer, and Herschel respectively, as well as near-UV 
(NUV) and FUV photometry observed primarily by 
GALEX and completed by the Ultraviolet /Optical Tele- 
scope (UVOT) on Swift. In the next sections, we de- 
scribe the observations and their reduction. The obser- 
vations were supplemented with mid-infrared (MIR) to 
FIR fluxes measured by the Infrared Astronomical Satel- 
lite {IRAS), the 70 /im and 160 /im Multiband Imag- 
ing Photometer (MIPS) detectors on Spitzer, and UBV 
fluxes from the Third Reference Catalog (RC3; de Vau- 
couleurs et al. 1991) where available in the literature 
through the NASA Extragalactic Database (NED). The 
MIPS 24 /im fluxes from these sources typically agree 
within the uncertainties with those we measure. Figures 
1-6 show from left to right the GALEX, 2MASS, Spitzer 
Infrared Array Camera (IRAC), and Herschel observa- 
tions of each galaxy. Some galaxies (e.g. NGC 3031 and 
M51A) have similar morphology from UV to FIR. In con- 
trast, others have distinct morphological differences be- 
tween the UV and IR, such as the FIR bright spots of 
NGC 2976 or the extended UV disk of NGC 3430. Ap- 
pendix C contains notes on the individual galaxies. 

3.1. Galaxy Distances 

All of the galaxies in our sample are nearby (within 60 
Mpc) and can therefore have peculiar velocities that con- 
tribute significantly to their recessional velocities. TuUy 
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TABLE 3 

Description of Herschel Observations 



Galaxy 


Instrument 


ObsID 




Date 


Obs. Mode 


Exposure (s) 




PACS Bands 




SPIRE 


1 S491 Q91 nfi 

J.O^Z J.yZ J.UU 


2010 


Mar 11 


-l-JcLig,e iVicL|J 


1076 












PACS 


1 S499n71 7n-7S 


2010 


Oct 26 


^r»nn A/Tot^ 
O'Neill iVicLJJ 


3624 


75 


/im, 


110 /im. 


170 /im 




OJr irvH/ 


1 Szl91 8'=%'=%S8 


zuuy 


wci o 


Large M!ap 










lo^zioDUoO-oO 


zuuy 


Or-f 1 7 


Scan M!ap 


99908 


75 


/xm, 


170 /im 




NGC3034 


SPIRE 




2009 


Oct 6 




2418 
















901 


IVrk^r 1 
iN ov xu 


Scan M!ap 


Ri=;zL9 

UO'IZ 


75 


/xm, 


170 /im 




NOC3n77 


SPIRE 


1342193015 


2010 


Mar 28 


Tj^rcp A/Tan 


2095 












PACS 




2011 


Mar 91 


^r*jin A/TciT* 


7356 


75 


)Lim, 


110 /im, 


170 /im 


NnCSI 8"^ /SI 87 /SI QD 


SPIRE 


1 S421 QfififiS 


2010 


Mav 1 S 


Tinycrfi A/InT* 
-Ljdigc iv±djj 


1035 










NnnSI R7 /SI QD 


PACS 


1 S499n71 


2010 


Oct 25 


OCcLll iVlcLjJ 


2708 


75 


)Ltm, 


110 /im, 


170 /im 


Nnr!S99fi/S997 


SPIRE 


1 S491 Q7S1 S 


2010 


Mpv so 


T.nTcr*::" A/TnT^ 
-LJdigC iVldjJ 


2624 












PAP^ 


1 SZL9991 1 zLfi ZL7 


901 
ZU lU 


iVldy -LO 


Scan M!ap 


1 1 7R 


75 


)Ltm, 


170 /im 






O-T i-Ttll/ 


1 '^zL99nQ9Rfi 


901 

ZU lU 


iN ov y 


Large M!ap 


QQQ 

yyy 












PACS 


1 S42221 1 04-07 


2011 


Mav 1 


>\r*JiTi A/Tfi'p* 
kjv_ydii iv±d|J 


2564 


75 


)Ltm, 


110 /im. 


170 /im 


l\TPP'^zl9zl /'^ZL'^n 


oJr irtii/ 


1 ^^491 Qt=;QzLfi 
lO^ZiyOy^D 


901 
ZU lU 


iviay o 


Large M!ap 


1 fil R 

ID lO 










MPP'^zl/LS /TTPPfifll 

iN 'or'^O^t^tO/ U'orV^OUiO 


QPTPF' 

O-T i-TLIl/ 


io4tz±ooooy 


900Q 

zuuy 


wet o 


Large Map 


1 5^*^^^ 
looo 










iN Vjv^ooyu/ iv^oy^t 


QPTPF' 

O-T i-tXil/ 


1 '^/L91 QQ'^/L/L 


901 
ZUiU 


Tim 9Q 

Jim zy 


Large Map 


4toy 












PA PQ 


iOfizzlUDUU-UO 


001 o 
zUiU 


Mi^Tr QO 

InOV oU 


Scan Map 


Z^tDz 


75 


)Lim, 


110 /im. 


170 /im 




PACS 


134211104-05 


2010 


Dec 13 


Scan Map 


524 


110 Aim 


, 170 /im 


NGC3786/3788 


SPIRE 


1342223233 


2011 


Jun 28 


Small Map 


169 










PACS 


1342223319-20 


2011 


Jun 29 


Scan Map 


104 


75 


)Lim, 


170 /im 




NGC4038/4039 


SPIRE 


1342188686 


2009 


Dec 29 


Large Map 


710 








PACS 


1342187836-39 


2009 


Dec 8 


Scan Map 


2662 


75 


)Lim, 


110 /im. 


170 /im 


NGC4618/4625 


SPIRE 


1342188755 


2009 


Dec 31 


Large Map 


1052 










NGC4625 


PACS 


1342210468-71 


2010 


Nov 19 


Scan Map 


2708 


75 


)Lim, 


110 /im. 


170 /im 


NGC4647/4649 


SPIRE 


1342188778 


2009 


Dec 31 


Large Map 


4295 










M51 


SPIRE 


1342188589 


2009 


Dec 26 


Large Map 


1577 












PACS 


1342188328-29 


2009 


Dec 20 


Scan Map 


4422 


75 


/im, 


170 /im 




NGC5394/5395 


SPIRE 


1342236140 


2012 


Jan 1 


Large Map 


1253 








PACS 


1342211285-88 


2010 


Dec 17 


Scan Map 


2200 


75 


/im. 


110 /im. 


170 /im 


MlOl 


SPIRE 


1342188750 


2009 


Dec 30 


Large Map 


9443 






PACS 


1342198471-74 


2010 


Jun 16-17 


Scan Map 


38077 


75 


/im. 


110 /im. 


170 /im 


NGC5474 


SPIRE 


1342188751 


2009 


Dec 26 


Large Scan 


1052 












PACS 


13422077178-81 


2010 


Oct 26 


Scan Map 


2708 


75 


/im. 


110 /im. 


170 /im 



TABLE 4 

Description of GALEX Observations 



Galaxy 


Tilename 


NUV 


FUV 






Date Exposure (s) 


Date Exposure (s) 



NGC2976 


GI2_024002_NGC2976_stream 


2006 


Jan 04 


18113, 


.55 


2006 


Jan 04 


17212.50 


NGC3031 


GI1_071001_M81 


2005 


Jan 12 


29421, 


.55 


2006 


Jan 05 


14706.70 


NGC3034 


NGA_M82 


2009 


Jan 31 


17311, 


.95 


2009 


Jan 31 


11527.35 


NGC3185/3187/3190 


NGA_NGC3190 


2004 


Jan 30 


3545 


.80 


2005 


Feb 19 


1299.15 


NGC3395/3396/3424/3430 


GI1_078004_NGC3395 


2006 


Mar 23 


2666 


.15 


2006 


Mar 23 


1500.10 


NGC3448/UGC6016 


AIS_92 


2004 


Feb 4 


423 


.00 


2004 


Apr 21 


143.00 


NGC3690/IC694 


AIS_99 


2007 


Feb 13 


211, 


.00 


2007 


Feb 13 


211.00 


NGC3786/3788 


AIS_111 


2007 


Feb 20 


103 


.05 


2007 


Feb 20 


103.05 


NGC4038/4039 


NGA_Antennae 


2004 


Feb 22 


1541, 


.30 


3004 


Feb 22 


1541.30 


NGC4618/4625 


NGA_NGC4625 


2004 


Apr 5 


3259 


.00 


2004 


Apr 5 


3259.00 


NGC4647/4649 


GI1_109003_NGC4660 


2005 


Apr 30 


3113 


.25 


2008 


Apr 23 


1624.10 


M51 


GI3_050006_NGC5194 


2007 


May 29 


10216 


.20 


2007 


May 29 


10216.20 


NGC5394/5395 


GIl_026018_Arp84 


2006 


Apr 12 


4268, 


.65 


2007 


May 30 


2811.40 


MlOl 


GI3_050008_NGC5457 


2008 


Apr 4 


13294 


.40 


2008 


Apr 4 


13293.4 


NGC5474 


NGA_NGC5474 


2003 


Jun 19 


1610 


.00 


2003 


Jun 19 


1610.10 
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TABLE 5 

Description of Swift UVOT Observations 









Exposure Times (s) 


Gcilcixy 


ObsID 


Date 


uvw2 


uvm2 


uvwl 


NGC3226/3227 


00031280001 


2008 Nov 4 


342 


249 


352 


NGC3226/3227 


00031280002 


2008 Nov 5 


704 


511 


346 


NGC3226/3227 


00031280003 


2008 Nov 12 


692 


424 


372 


NGC3226/3227 


00031280004 


2008 Nov 13 


744 


538 


372 


NGC3226/3227 


00031280005 


2008 Nov 21 


744 


522 


381 


NGC3226/3227 


00031280006 


2008 Nov 22 


763 


137 


381 


NGC3226/3227 


00031280007 


2008 Nov 25 


763 


531 


381 


NGC3226/3227 


00031280008 


2008 Nov 27 


763 


196 


246 


NGC3226/3227 


00031280009 


2008 Dec 2 








293 


NGC3226/3227 


00031280010 


2008 Dec 3 


274 


349 


126 


NGC3424/3430 


00091132001 


2011 Apr 16 








1976 


NGC3424/3430 


00091132003 


2011 Jun 28 











NGC3424/3430 


00091132004 


2011 Jul 4 





80 


1315 


NGC3424/3430 


00091132005 


2011 Jul 7 


302 








NGC3424/3430 


00091132006 


2011 Jul 8 





1877 





NGC3424/3430 


00091132007 


2011 Oct 7 


750 


988 





NGC3424/3430 


00091132008 


2011 Oct 10 












Note. — The Swift observation ID number (Col. 2) and the start date of 
each observation (Col. 3) are given for each observation of each object for 
which observations with minimal coincidence losses exist. Exposure times in 
the each filter are given in Col. 4-6. 



NGC297G 




Fig. 1. NGC 2976, NGC 3031, NGC 3034, and NGC 3077 (from top to bottom) as observed, from left to right, by GALEX (NUV in 
yellow; FUV in blue), 2MASS (J in blue, H in green, and Ks in red), IRAC (3.6 jim in blue, 4.5 fim in green, and 8.0 jim in red), and 
Herschel (FAGS 75 /im in blue, FAGS 170 /im in green, and SFIRE 250 /im in red). The longer wavelength IRAG observations of NGG 
3034 were saturated, so 4.5 yum is shown in yellow instead. NGG 3077 was not observed by either GALEX or Swift. At the distance of 
these galaxies, l';^l.l kpc. 
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Fig. 3. As Figure 1, but for NGC 3690/IC694, NGC 3786 (bottom)/NGC 3788 (top), and NGC 4038/4039. The 8 fim IRAC image of 
NGC 3690/IC 694 is saturated in the nuclei of the two galaxies, resulting in the blue-green artifacts. At the distance of these galaxies, 1' 
is approximately 14 kpc (NGC 3690), 11-12 kpc (NGC 3786/3788), and 7.4 kpc (NGC 4038/4039). 
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NGC3395/9G 



..m 




-a ' 




Fig. 4.— As Figure 1, but for NGC 3226 (upper)/NGC 3227 (lower), NGC 3395 (right)/NGC 3396 (left), NGC 3424 (riglit)/NGC 3430 
(left), and NGC 3448 (left)/UGC 6016 (right). NGC 3226/3227 was not observed by GALEX but by Swift Their left image show the Swift 
observations through the UVWl filter in blue, the UVM2 filter in green, and the UVW2 filter in red. NGC 3424/30 and NGC 3448/UGC 
6016 were not observed with FAGS, so the right image only shows the SFIRE 250 fim image as in Figure 2. UGC 6016, while having 
significant extended diffuse emission in the UV, is not well detected in the IR bands. At the distance of these galaxies 1' is approximately 
6-8 kpc. 



et al. (2008) recently compiled redshift-independent dis- 
tances for nearby galaxies with velocities less than 3000 
km s~^ using alternate methods including Cepheids 
(Freedman et al. 2001), the luminosity of stars at the 
tip of the red giant branch (Karachentsev et al. 2006), 
surface brightness fluctuations (Tonry et al. 2001), and 
the Tully-Fisher relation (Tully & Fisher 1977). Dis- 
tances to additional galaxies based on their group or clus- 
ter association are given in the Extra-galactic Distance 
Database^^ (EDD; Tully, 2010, private communication). 

http://edd.ifa.hawaii.edu 



Twenty-six of our galaxies have distance moduli given by 
either Tully et al. (2008), Tully (1994), or EDD. For the 
five galaxies lacking distance moduli, we obtained helio- 
centric velocities from the PSCz catalog (Saunders et al. 
2000; NGC 3690/IC 694, NGC 5394, and NGC 5395) 
and RC3 (UGC 6016), which we corrected to account 
for the velocity field of Virgo, the Great Attractor, and 
the Shapley supercluster, following Mould et al. (2000). 
Distances were then calculated assuming Hq =72 km s~^ 
Mpc~^. The distances are given in Table 1. 
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Fig. 5. As Figure 1, but for NGC 4618, NGC 4625, NGC 4647 (right)/NGC 4649 (left), and M51. The right images of NGC 4618 and 
NGC 4647/49 only show the SPIRE 250 /xm image. At the distance of these galaxies, 1' is approximately 2.1-2.4 kpc (NGC 4618/4625), 5 
kpc (NGC 4647/4649), and 2.2 kpc (M51). 



3.2. Infrared Photometry 
3.2.1. Spitzer Observations 

The IRAC (Fazio et al 2004) and MIPS (Rieke et al. 
2004) 24 jam observations were taken as part of a vari- 
ety of programs, including the main SIGS program (PID 
20140; P.I. A. Zezas), which also observed galaxy groups 
that had not previously been observed. The observa- 
tion parameters are given in Table 2. The IRAC Basic 
Calibrated Data (BCD) were retrieved from the Spitzer 
archive and cleaned before being coadded into mosaics 
with 0.6^' pixels using IRACproc (Schuster et al. 2006). 
The MIPS 24 jam BCDs were merged to form mosaics 
with 2.45'^ pixels using the Mosaicker and Point Source 
Extractor package (MOPEX; Makovoz & Khan 2005). 



The reduction of these data will be described in detail 
in Brassington et al. (2013, in preparation). While the 
pipeline versions range from S13-S18, the difference be- 
tween the pipelines are minor and do not impact signif- 
icantly the photometry^^. The pipeline version for each 
galaxy is given in Table 2. 

3.2.2. Herschel Observations 

The parameters for the Herschel SPIRE (Griffin et al. 
2010) and PACS (Poglitsch et al. 2010) observations are 
given in Table 3. The Herschel data were taken as part 
of two Science Demonstration Phase programs (P.I.s C. 

■^'^ http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/ 
iracinstrument handbook /79/ 
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Fig. 6. As Figure 1, but for NGC 5394 (top)/NGC 5395 (bottom), MlOl, and NGC 5474. At the distances of these galaxies, 1' is 
approximately 16 kpc (NGC 5394/5395) and 1.7-2.0 kpc (MlOl/NGC 5474). 

Wilson and S. Eales), four Key Project programs (P.I.s 
R. Kennicutt, S. Eales, C. Wilson, and E. Sturm), and 
one Guaranteed Time program (P.L L. Spinoglio). All of 
the galaxies were observed by SPIRE at 250 /im, 350 /im, 
and 500 /im; this was part of the selection criteria of this 
sample. Approximately 50% of the sample were observed 
in all three PACS bands and an additional '^25% were 
observed at 75 /im and 170 /im. 

The data were retrieved from the Herschel Science 
Archive and processed using the calibration trees of ver- 
sion 8.0.1 of the Herschel Interactive Processing Envi- 
ronment (HIPE; Ott 2010). This processing was ac- 
complished using the default pipeline scripts available 
through HIPE to make Large Map mode mosaics for the 
SPIRE data and extended source mosaics with MadMap 
for PACS data. We discuss additional details regarding 
the processing of PACS data in Appendix Bl. 



3.2.3. 2MASS Observations 

NIR mosaics of the sample galaxies observed as part of 
the 2MASS (Skrutskie et al. 2006) were retrieved from 
the NASA/IPAC Infrared Science Archive^^, and from 
the Large Galaxy Atlas (Jarrett et al. 2003) when possi- 
ble. The counts measured in the images were converted 

NASA/IPAC Infrared Science Archive is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under 
contract with NASA. 



to Janskys using the zero points of Cohen et al. (2003). 
We compared our fluxes measured in the apertures de- 
scribed in §3.4 to the total fluxes given in NED from 
Jarrett et al. (2003) and the 2MASS Extended Object 
Catalog and found good agreement. 

3.2.4. Ancillary IRAS Photometry 

IRAS photometry was obtained from the HIRES At- 
las (Surace et al. 2004), the IRAS Revised Bright Galaxy 
Sample (Sanders et al. 2003), the IRAS Bright Galaxy 
Sample (Soifer et al. 1989), and the Faint Source Cat- 
alogue (Moshir et al. 1990). The latter three catalogs 
present photometry derived from the native IRAS beam 
size of 2^ — 5'; this can be problematic for systems in close 
interaction phases. We therefore preferentially used the 
HIRES Atlas, which was reprocessed with 30'^ — 1.^5. In 
the one system where only low-resolution photometry is 
available and the galaxies are close enough for contami- 
nation to occur, we do not include the IRAS photometry 
in our analysis. 

3.3. Ultraviolet Photometry 

3.3.1. GALEX Observations 

Twenty-eight of our sample galaxies were observed by 
GALEX] three sources within the sample (NGC 3226, 
NGC 3227, and NGC 3077), however, were not observed 
due to the presence of nearby bright stars. For the galax- 
ies with GALEX photometry, mosaics of the longest ob- 



12 



Lanz et al. 



servations were retrieved from the Mikulski Archive for 
Space Telescopes^^ using GalexView version 1.4.6. The 
details of those observations are given in Table 4. The 
NUV observation of NGC 3690/IC 694 was reprocessed 
by D. Neill at our request to correct a masking prob- 
lem. We use the conversions from count rate to fluxes 
provided by Goddard Space Flight Center (2004)^^. 

3.3.2. Swift UVOT Observations 

Most of the gaps in the GALEX coverage can be filled 
in with data from the Swift UVOT telescope, which has 
three UV filters that bracket the GALEX NUV filter 
in mean wavelength. Two of the three galaxies lacking 
GALEX data, NGC 3226 and NGC 3227, were observed 
by UVOT. Unfortunately, NGC 3077's nearby bright star 
exceeded the tolerances of this telescope as well. We orig- 
inally planned to use existing UVOT photometry for all 
our sample. We obtained the raw data and exposure 
maps from the Swift archive for the seventeen galax- 
ies with UVOT data and coadded the observations into 
one mosaic and exposure map per UV filter per interact- 
ing system. However, as described by Hoversten et al. 
(2011), the photon-counting nature of the Swift detec- 
tors makes them vulnerable to coincidence losses, which 
become significant when the count rate is greater than 
0.007 counts per second per pixel. We calculated count 
rate maps to determine where coincidence losses need to 
be taken into account. Due the difficulties associated 
with coincidence losses in extended sources, described in 
greater length in Appendix B2, we opted only to use the 
UVOT data for the missing GALEX objects NGC 3226 
and NGC 3227. We added one test case, NGC 3424, to 
confirm that the UVOT data yielded fiuxes consistent 
with GALEX and found good agreement. The details 
of the observations of these three galaxies are given in 
Table 5. To convert the count rate to fluxes, we used 
the conversion assuming a stellar spectrum described in 
Breeveld et al (2010). 

3.4. Aperture and Uncertainty Determination 

For consistency, we used matched apertures across all 
wavebands in our photometric analysis. Generally, the 
IR emission of galaxies is more extended that their UV 
emission. However, some of the galaxies are more ex- 
tended in the UV than in the IR (e.g. NGC 3430). We 
used the SExtractor algorithm (Bertin & Arnouts 1996) 
to determine Kron apertures in both the NUV and the 
3.6 /xm IRAC images. In all cases, the larger of the two 
apertures was then used to measure the integrated galaxy 
flux at all wavelengths in order to obtain flux from a con- 
sistent area of each galaxy across our wavelength range. 
The size and position angle of each aperture as well as 
on which image it was determined is given in Table 1. 
Background regions were selected to mimic the content 
of background and foreground objects in the aperture 
on the outskirts of the galaxies. Once the aperture was 
selected, flux densities in the aperture and background 

STScI is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS5-26555. 
Support for MAST for non-HST data is provided by the NASA Of- 
fice of Space Science via grant NNX09AF08G and by other grants 
and contracts. 

http://galexgi.gsfc.nasa.gov/docs/galex/FAQ/counts_backgro 



regions were measured using the analysis tools of the 
SAOImage DS9 (Joye & Mandel 2003). 

Due to the proximity of some members of the same in- 
teracting system, their apertures can overlap. We dealt 
with these situations in one of three ways. For sig- 
nificantly overlapping systems (NGC 4038/4039, NGC 
3690/IC 694, and NGC 3395/3396), separate apertures 
could not be robustly determined. In these cases, we 
treated the combined system as a single object. Second, 
there were two systems (M51 A/B and NGC 5394/5495) 
where the aperture for the smaller galaxy was mostly 
contained within the aperture of the larger galaxy, but 
it was clear that the emission in the overlap area came 
from the smaller galaxy. In these cases, we subtracted the 
emission and area of the overlap region from that of the 
larger aperture. Third, there were three systems (NGC 
3226/3227, NGC 3786/3788, and NGC 4647/4649) where 
the aperture overlapped but without significant contami- 
nation. In these cases, we extrapolated the expected flux 
in the overlap area from the surface brightness in the rest 
of the elliptical aperture at the same radii. 

The Spitzer fiuxes required aperture corrections. We 
determined the effective radius of the elliptical aper- 
ture^^ and used the extended source flux corrections 
given in the IRAC Instrument Handbook^^. For the 
MIPS 24 /im aperture corrections, we interpolated be- 
tween the aperture corrections given in the MIPS Instru- 
ment Handbook^^. The GALEX data were corrected for 
obscuration due to Milky Way dust using the extinction 
laws given by Wyder et al (2005). 

Uncertainties in the absolute fluxes are the sum in 
quadrature of a statistical uncertainty and a calibration 
uncertainty. The Spitzer bandpass uncertainties are typ- 
ically dominated by the calibration uncertainty of 3% for 
IRAC (Cohen et al 2003) and 4% for MIPS 24 /im (En- 
gelbracht et al 2007). We used a calibration uncertainty 
of 10% for the GALEX data (Goddard Space Flight Cen- 
ter 2004) and a 5-15% uncertainty for the Swift bands 
(Poole et al 2008), and the statistical uncertainty is cal- 
culated using Poisson statistics. We used a 7% calibra- 
tion uncertainty for the SPIRE bandpasses (Swinyard et 
al 2010) and 10% for the PACS bandpasses (Paladini et 
al 2012) and followed Dale et al (2012) in calculating 
the statistical uncertainty. The photometry results for 
GALEX, Swift, and 2MASS; Spitzer; and Herschel are 
provided in Table 6-8, respectively. When flux is not de- 
termined significantly, we provide 3(J upper limits, but 
we do not provide lower limits in cases of saturated im- 
ages. The additional photometry from the literature is 
given in Table 9. 

4. SED FITTING WITH MAGPHYS 

4.1. Fitting Process 

To estimate SFR, specific star formation rates (sSFRs) , 
stellar and dust masses, and dust temperatures, we used 
the SED fitting code MAGPHYS (da Cunha et al 2008). 
MAGPHYS fits SEDs with a combination of UV-NIR 
stellar spectral libraries from Bruzual & Chariot (2003) 

leff = ^/ab for semi-major axis a and semi-minor axis h 
http:/ /irsa.ipac. caltech.edu/data/SPITZER/docs/irac/ 
iracinstrument handbook/30/ 

http:/ /irsa.ipac. caltech.edu/data/SPITZER/docs/mips/ 
, .hti]ii]bsinstrumenthandbook/50/ 
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Fig. 7. — SEDs for NGC 2976 (top) and NGC 3031 (bottom) with data shown as red points, the best fit model plotted in black, and the 
stellar emission in the absence of dust shown in blue. The components of the infrared emission are over plotted: PAH emission (dotted, 
green line), MIR emission at 130 K and 250 K (red, dashed line), warm 30-60 K dust emission (dot-dashed, purple line), and cold 15-25 
K dust emission (long dashed, orange line). Below the fitted SEE) is plotted the fractional difference between the model and data. To the 
right of the SED, we plot a subset of the probability distribution functions (PDFs) of the fitted parameters for (from left to right): stellar 
mass, dust luminosity, and dust mass (top) and SFR, sSFR, and cold dust temperature (bottom). 



and a simple, physically- motivated model for IR emis- 
sion from dust developed in da Cunha et al. (2008). It 
models the ISM as a mix of diffuse dust interspersed 
with denser, warmer stellar birth clouds. MAGPHYS 
also includes a set of UV— NIR libraries that modify 
the Bruzual & Chariot (2003) population synthesis with 
the Bruzual (2007) population synthesis, which provides 
different treatment of post-Asymptotic Giant Branch 
(AGB). We fit our UV to FIR SEDs with MAGPHYS 
with and without the post- AGB modifications and found 
consistent results; from here on, we only use the results 
with the earlier Bruzual & Chariot (2003) libraries as 
their treatment of the post- AGB stars is more consistent 
with current understanding (e.g., Zibetti et al. 2012). 
The IR dust libraries have five components: a fixed poly- 
cyclic aromatic hydrocarbon (PAH) spectrum shape de- 
rived from the Ml 7 SW star- forming region (Madden et 
al. 2006), a NIR continuum associated with the PAH 
emission modeled by a modified blackbody {j3 = 1) at 
850 K, a hot MIR continuum modeled by the sum of two 
modified blackbodies {/3 = 1) at 130 K and 250 K, a 
warm (30-60 K) dust component modeled as a modified 
blackbody (/3 = 1.5), and a cold (15-25 K) dust com- 



ponent modeled as a modified blackbody {P = 2). The 
warm dust component is assumed to exist both in the dif- 
fuse ISM and in denser birth clouds, while the cold dust 
exists only in the diffuse ISM. MAGPHYS determines 
probability distribution functions (PDFs) for the fitted 
parameters by combining UV— NIR and IR spectral li- 
braries such that the energy absorbed in the UV/ visible 
regime is re-emitted in the IR. It gives both the best- 
fit obscured SED and the associated unobscured stellar 
SED. 

We input the measured and literature fluxes in our set 
of 28 filters to MAGPHYS and examined the variation 
in derived galaxy properties including dust luminosity, 
SFR, sSFR, stellar and dust mass, and dust tempera- 
tures and discuss the results below. Then we performed 
five additional fits: one without UV photometry, one 
without SPIRE photometry, one without any photom- 
etry at wavelengths A > 30 /im, one without either UV 
or SPIRE photometry, and one with only UBV, 2MASS, 
and IRAC photometry. For these fits for each galaxy, 
we examined the median and 68% confidence interval for 
the fitted parameters to determine the influence of the 
particular dataset on the value of (and constraints on) 
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Fig. 8. As Figure 7, but for NGC 3034 (top), NGC 3077 (middle), and NGC 3185 (bottom). 



these properties, as we elaborate below. For the non- 
interacting galaxies, we only performed the fit with all of 
the available photometry. 

4.2. SED Fits 

Figures 7—16 show the SEDs along with the best-fit 
models for our interacting galaxies. The contributions of 
the different components of the IR model described above 
are also shown. The median and 68% confidence interval 
of these parameters are given in Table 10. Just as the 



galaxies exhibit a variety of UV versus IR morphologies 
(Figures 1—6), the SEDs have a corresponding range of 
relative UV, NIR, and FIR emission. For example, some 
galaxies (e.g., NGC 3190 or M51B) have very little UV 
flux in comparison with their visible and IR flux, while 
others (e.g., NGC 3187) have almost as much UV flux 
as IR flux. The SEDs also show a range of obscuration 
from the heavily obscured galaxies (e.g., NGC 3690) to 
relatively unobscured galaxies (e.g., NGC 4618), or ones 
with extended UV disks such as NGC 3430 or UGC 6016. 
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Fig. 9. As Figure 7, but for NGC 3187 (top), NGC 3190 (middle), and NGC 3226 (bottom). 



Appendix C briefly describes each galaxy, discussing any 
particular issues regarding the photometry and the SEE) 
fitting. Note in particular that fits to edge-on galaxies 
tend to over-estimate the amount of UV obscuration and 
hence the model UV fluxes tend to be too low compared 
to the observations. 

5. DISCUSSION 
5.1. Variation in SED Shape with Interaction Stage 



We now discuss the shape of the SED as a function 
of the interaction stage. We first normalized each SED 
to its 2MASS Ks luminosity. Emission in the 2MASS 
Ks filter is dominated by the old stellar populations and 
hence is a good proxy for stellar mass. The comparison 
between the SEDs is shown in Figure 17. In the inset, we 
show the median normalized SED for each class of galax- 
ies. The SED shapes between the three stages vary by 
approximately as much as the variations within a stage. 
However, there are some significant variations, especially 
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As Figure 7, but for NGC 3227 (top), NGC 3395/3396 (middle), and NGC 3424 (bottom). 



Fig. 10. 

in the Stage 4 SEDs compared to the Stage 2 and 3 SEDs. 
Stage 4 galaxies typicahy have more emission from the 
hot/warm dust than earher interaction stages, as evi- 
denced by the stronger 10—60 /im emission relative to 
their stellar mass. Further, Stage 4 galaxies tend to have 
more warm dust relative to their cold dust FIR emission. 
In contrast to this variation in the relative MIR emission, 
all three stages have similar ratios of NIR stellar emis- 
sion to FIR emission from cold dust. We will discuss the 
statistical significance of these results in the next section. 



These differences in the SEDs are consistent with the 
results of simulations, which predict that an integrated 
SED of an interacting system becomes hotter at merger 
coalescence during the peak of starburst and AGN activ- 
ity (e.g., Hayward et al. 2011, 2012a; Younger et al 2009; 
Narayanan et al. 2010a, 2010b). However, the increase 
in temperature in this sample of galaxies is unlikely to 
be driven purely by AGN activity. None of the galaxies 
in our sample have typical AGN colors, based on their 
location outside the region in IRAC color space defined 
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Fig. 11. As Figure 7, but for NGC 3430 (top), NGC 3448 (middle), and UGC 6016 (bottom). Note that the axes of the UGC 6016 
plots have smaller values than the rest of the plots. 



by Stern et al (2005) (Figure 18). Further, while five of 
the sample galaxies are classified as Seyfert galaxies and 
three are classified as low-ionization nuclear emission- 
line region (LINER) galaxies, they are found in all three 
stages. We used the software DECOMPIR^^ (Muhaney 
et al. 2011) to estimate the AGN contribution to the 
8—1000 jiva and 8—35 /im emission for these nine galax- 

http:sites.google.com/site/decompir 



ies based on the 8-500 /im photometry. We give the 
individual contributions in the descriptions in the Ap- 
pendix C; the range of the contribution to the total IR 
is up to 10% with some larger, more uncertain values up 
to 25%. Further, we do not find significant differences 
in the AGN contribution to either the total infrared lu- 
minosity or mid-infrared luminosity, where the Stage 4 
SEDs are typically brighter, between the stages. Hence, 
the effect of AGN on the SEDs of the sample galaxies is 



22 



Lanz et al. 




1.0 



1 0.0 1 00.0 

(obaerved -frame] 



1 000 .0 



1.0 



1 0.0 1 00.0 

(observed -frame] 



100Q.0 



1 .0 1 0.0 1 00.0 

[observed -frame) 



IQOO.Q 



g 10 11 



ID 11 12 5 6 7 a 



-7 2 4-12-11-10-9 



la ie 20 32 24 




9 10 11 



9 ID 11 



5 6 7 a 



-7 2 4-l3-l?-1 1-10 



16 ie 20 32 24 




9 10 11 



9 ID 11 



5 6 7 a 



-7 2 



-13-12-11-10 



16 16 20 22 24 



Fig. 12. As Figure 7, but for NGC 3690/IC 694 (top), NGC 3786 (middle), and NGC 3788 (bottom). Note that the axes of the NGC 
3690 plots have larger values than the rest of the plots. 

modest and does not affect our conclusions. 

The more strongly interacting systems demonstrate a 
tendency to have, on average, younger stellar popula- 
tions, resulting in stronger UV emission relative to their 
NIR emission. Since the UV emission has only been cor- 
rected for Milky Way extinction, additional intrinsic ex- 
tinction could increase this effect. Ignoring NGC 4649, a 
large elliptical that has very little MIR-FIR emission, the 
UV bands reflect this tendency and show a large amount 



of variation between galaxies, which is likely due to the 
different amounts of dust attenuation and the sensitivity 
of the UV to recent star formation history. Comparing 
to the stellar mass proxy of 2MASS Ks^ Stage 4 galaxies 
typically have a 1:2 luminosity ratio between emission in 
the GALEX bands and 2MASS Ks band, whereas Stage 
2 galaxies typically have a 1:10 luminosity ratio. 
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13. As Figure 7, but for NGC 4038/4039 (top), NGC 4618 (middle), and NGC 4625 (bottom). 



Fig. 

5.2. Variations in Galaxy Parameters with Interaction 

Stage 

Figures 19-20 show the histograms and cumulative dis- 
tributions for the three interaction stages defined in §2.2 
for four parameters whose PDFs are shown in Figures 
7-16 compared to the non-interacting galaxies. Table 11 
gives the results of the Kolmogorov-Smirnov (K-S; e.g., 
Press et al 2007) tests performed to determine the like- 
lihood that the samples for the different stages originate 



from a common parent population for each parameter. 
With the number of galaxies in our sample, we do not 
have the same statistical power that the full SIGS sample 
will have (Brassington et al. 2013, in prep.). 

The mass of dust is one of four parameters derived from 
the SED that differs with marginal statistical significance 
between the Stage 3 and Stage 4 galaxies {p = 0.02). 
Both the mass of warm dust and of cold dust likewise 
differ. The warm dust mass also differs between the 
Stage 4 galaxies and the non-interacting galaxies. In 
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Fig. 14. As Figure 7, but for NGC 4647 (top), NGC 4649 (middle), and M51A (bottom). 



contrast, the distribution of dust mass in the ISM of 
non-interacting galaxies fahs between the distributions 
of Stage 2/3 and Stage 4 galaxies. The dust mass in the 
ISM and stellar mass also each show only small increases 
in their median values as the interaction sequence pro- 
gresses, and their cumulative distributions are not signif- 
icantly different. Indeed, both the stellar mass {p = 0.96) 
and total dust mass {p = 0.96) for Stage 2 and Stage 3 
galaxies are consistent with coming from the same popu- 
lations. These trends are consistent with the SED shapes 



described in §5.1, where the Stage 4 SEDs typically show 
more emission from the hot /warm dust relative to the 
cold dust emission and stellar emission. 

The dust luminosity (Figure 19, left) shows marginally 
statistically significant differences between Stage 3 and 
Stage 4 galaxies {p = 0.01) and between the non- 
interacting galaxies and the Stage 4 galaxies {p = 0.02). 
The median dust luminosity increases with interaction 
stage by over an order of magnitude between Stage 2 
and Stage 4. Elbaz et al. (2011) defined an "IR main 
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15. As Figure 7, but for M51B (top), NGC 5394 (middle) and NGC 5395 (bottom). 



Fig. 

sequence" of galaxies in which the ratio of total IR lumi- 
nosity to 8 /im luminosity has a Gaussian distribution. 
We examined this ratio for our sample and found good 
agreement with the expected distribution. The only in- 
teracting galaxy that lies off this relation, by roughly 
an order of magnitude, is the large elliptical NGC 4649, 
as would be expected. Elbaz et al. (2011) also defined 
two modes of star formation: a "normal" mode exem- 
plified by the galaxies on the IR main sequence and a 
"starburst" mode with excess sSFR in comparison. Our 



sample's agreement with the IR main sequence indicates 
that our set of interacting galaxies do not contain sys- 
tems with significantly increased sSFR. 

We also considered the evolution of the cold and warm 
dust temperatures. The cold dust temperature (Figure 
19, right) is the third parameter showing evidence for 
differences between the Stage 3 and Stage 4 {p = 0.01) 
and between the non-interacting galaxies and the Stage 4 
galaxies {p = 0.01). The cold dust temperature's median 
value varies in a similar manner to the IR luminosity, in- 
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Fig. 16. As Figure 7, but for MlOl (top) and NGC 5474 (bottom). 

creasing between Stages 3 and 4 but relatively constant 
between Stages 2 and 3; the median value of the warm 
dust temperature is by contrast fairly constant. Only 
Stage 4 does not span the range of the 15—25 K cold dust 
temperature, while in Stage 2 and Stage 4 the warm dust 
temperatures are confined to the 45—60 K range. The 
similarity in evolutionary trend in the IR luminosity and 
cold dust temperature is likely due to the predominance 
of cold dust mass and luminosity in the total dust esti- 
mates. We might expect a similar trend to be exhibited 
in the temperature of the warm dust primarily present in 
the stellar birth clouds, however the warm dust tempera- 
ture is less well constrained in MAGPHYS than the cold 
dust temperature. The cold dust contribution and warm 
dust temperature are both correlated with the warm dust 
mass (which drives the MIR continuum intensity and the 
shape of the SED). Therefore, the warm dust tempera- 
ture PDFs tends to be broader, with a 68% range that is 
typically four to five times the size of the cold dust 68% 
confidence interval. However, since the cold dust mass 
is typically over ~ 80% of the total dust mass, the total 
dust mass is still fairly well constrained. 

The SFR, shown in Figure 20 (left), shows an increase 
in median value with interaction stage, like the the dust 
mass and dust luminosity, and has a marginally signif- 
icant probability the same population did not yield the 
Stage 4 galaxies as well as the Stage 2 (p = 0.03), Stage 
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3 {p = 0.02), and non-interacting {p = 0.02) galaxies. 
Since the warm (30—60 K) dust, primarily heated by 
young stars with ages less than 10 Myr, contributes sig- 
nificantly to the total dust luminosity, it makes sense that 
these three parameters show similar evolution over the 
interaction stages. However, an increase in SFR between 
the different stages could be also attributed to our Stage 

4 galaxies simply being larger with greater gas reservoirs. 
To test this, we also examined the evolution of sSFR over 
the interaction stages (Figure 20, right). In contrast to 
SFR, we do not find much difference in the median values 
of the sSFR, and the cumulative distributions are very 
similar in both width and normalization. We therefore do 
not see enhanced sSFR in more evolved mergers, consis- 
tent with the ratios of total IR luminosity to 8 /am lumi- 
nosity being due to a "normal" mode of star formation for 
our whole sample. We also do not see differences between 
the distributions of the sSFR of the three interacting 
galaxy samples and of the non-interacting galaxies. Xu 
et al. (2010) and Yuan et al. (2012) examined star forma- 
tion in a sample of local major mergers and both found 
that that the sSFR distributions of galaxies in spiral- 
spiral interactions and non-interacting systems were un- 
likely to originate from the same population (based on a 
K-S test: p = 0.03 — 0.04). However, they also found a 
mass dependence in the enhancement of sSFR; only those 
galaxies with > were found to have signif- 
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TABLE 10 
MAGPHYS Parameters 
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Note. — Median of the PDFs determined by MAGPHYS. The uncertainty given is the 68% confidence range. 
Averaged over 100 Myr 



icant enhancements. Our sample shows a similar trend, 
but we only have six galaxies with > lO^^-^M© of 
which one is an elliptical, and they are spread across 
Stages 3 and 4. 

We also examined the star formation efficiency, which 
we define as the ratio of SFR to warm (30—60 K) birth 
cloud dust mass, a proxy for molecular gas mass. While 
the ratio of dust mass to gas mass is not necessarily the 
same between galaxies, this ratio provides a means of es- 
timating the star formation efficiency under the assump- 
tion of constant gas-to-dust ratio. This star formation 



efficiency allows us to compare the SFR taking into ac- 
count the variable gas reservoirs. We find that regard- 
less of stage, the star formation efficiency ranges over 
more than three orders of magnitude and the cumulative 
distributions show no evidence of originating from dif- 
ferent populations. This result agrees with the findings 
of Casasola et al. (2004) who found similar star forma- 
tion efficiencies, defined as the ratio of FIR luminosity (a 
proxy for SFR) and molecular hydrogen mass, for inter- 
acting and non-interacting galaxies. 
All the apparent variations with interaction stage come 
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TABLE 11 
K-S Probabilities 



Parameter 


With Non-Interacting 
Stage 2 Stage 3 Stage 4 


Stages 2 and 3 


Between Stages 
Stages 3 and 4 


Stages 2 and 4 


Dust Luminosity 


0.238 


0.224 


0.017 


KJ. ID 1 


U.UUo 


n 1 on 


T^z-kfol T^nof A/Tooo 

loiai uusi iviass 


0.095 


0.072 


0.462 


u.yo4t 






Stellar Mass 


0.366 


0.655 


0.569 


0.964 


0.271 


0.129 


SFR 


0.980 


0.983 


0.017 


0.964 


0.022 


0.028 


sSFR 


0.569 


0.953 


0.606 


0.964 


0.768 


0.883 


Cold Dust Temperature 


0.337 


0.417 


0.006 


0.767 


0.008 


0.423 


Warm Dust Temperature 


0.569 


0.237 


0.106 


0.271 


0.271 


0.129 


Cold Dust Mass 


0.095 


0.072 


0.462 


0.964 


0.022 


0.423 


Total Warm Dust Mass 


0.682 


0.678 


0.019 


0.767 


0.022 


0.129 


Birth Cloud Warm Dust Mass 


0.912 


0.678 


0.019 


0.492 


0.022 


0.129 


SFR/Dust Luminosity 


0.238 


0.224 


0.999 


0.492 


0.271 


0.883 


SFR/Birth Cloud Warm Dust Mass 


0.119 


0.478 


0.462 


0.492 


0.492 


0.883 


Dust Mass/Stellar Mass 


0.644 


0.916 


0.857 


0.767 


0.271 


0.423 


Birth Cloud /Cold ISM Dust Mass 


0.719 


0.224 


0.037 


0.767 


0.271 


0.129 



Note. — The probabilities given in this table are the probability that the same population yielded the compared samples, so 
small values indicate a common population is unlikely. The marginally significant differences are given in bold. The total dust 
mass is composed of warm dust in the stellar birth clouds and the diffuse ISM and cold dust in the diffuse ISM. 




Fig. 17. — SEDs normalized to the 2MASS Ks band luminosity 
for each galaxy. Stage 2 (weakly interacting), Stage 3 (moderately 
interacting). Stage 4 (strongly interacting), and non-interacting 
galaxies are shown respectively in blue squares joined with solid 
lines, green triangles joined with long dashes, red diamonds joined 
with short dashes, and black circles joined with dotted lines, respec- 
tively. Inset, we show the median SED for each class of galaxies. 
These SEDs show a tendency for Stage 4 galaxies to have more 
hot-warm dust emission in the 10-60 /im range relative to both its 
cold dust emission in the SPIRE bandpasses and its stellar NIR 
emission, whereas the ratio of NIR to FIR emission is relatively 
consistent. Additionally, the more strongly interacting galaxies 
typically have a younger stellar population than the Stage 2 galax- 
ies as suggested by the relative amounts of UV to NIR emission. 
The two labeled galaxies are elliptical galaxies. 

with a few caveats. First, with only thirty-one galaxies, 
our sample has limited statistical power to identify signif- 
icant variations, especially with half the galaxies in Stage 
3. Analysis of trends in SFR and sSFR with the fuh SIGS 
sample (Brassington et al. 2013, in preparation), which 
covers the stages much more uniformly, will have greater 
statistical power (albeit these comparisons lack the Her- 
schel SPIRE data and SED analysis that provides more 
accurate SFR measurements). Second, our classification 
scheme permits us to examine parameter variations with 
respect to the strength of the interaction. While this se- 
quence crudely mimics an interaction, the dynamics of 
two colliding galaxies often includes multiple encounters 
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Fig. 18. — IRAC color-color plot of the SIGS galaxies (square) 
with our sample galaxies (diamonds) over-plotted, showing that 
none fall in the wedge identified by Stern et al. (2005) as galaxies 
hosting AGN. Our sample galaxies identified optically as Seyferts 
or LINERs in Keel et al. (1985) or Ho et al. (1997, 2000) and are 
marked with red and blue stars, respectively. 

prior to final coalescence, modifying the level of star for- 
mation at intermediate stages (e.g., Torrey et al. 2012), 
as well as the intensity of the final burst (e.g., Hopkins 
et al. 2008, 2009). As a result, interacting systems often 
do not progress linearly through the interaction stages 
defined by our classification system. 

5.3. Relative Importance of Specific Data Sets in 
Constraining Galaxy Parameters 

For each galaxy, we ran six MAGPHYS fits to measure 
the relative importance of UV, SPIRE, and MIR-FIR 
data in constraining the value of the derived SFR; stel- 
lar mass; sSFR; and dust temperatures, luminosity, and 
masses. We did this by comparing fits with all available 
data with fits using a subset of the complete dataset in 
order to determine if and how the absence of a particular 
dataset results in a systematic over- or under-estimation 
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Fig. 19. — Histograms (top) and cumulative distributions (bottom) of the dust luminosity (left) and cold dust temperature (right) derived 
by MAGPHYS for the three interaction stages defined by the classification system of Dopita et al. (2002), where Stages 2-4 are weakly 
(blue, solid), moderately (green, long dashed), and strongly (red, short dashed) interacting, respectively. Non-interacting galaxies are 
shown in black dot-dashed lines. There is a difference in both luminosity and temperature between the Stages 3 and 4 populations and 
between Stage 4 and non-interacting populations as defined by a K-S Test (see Table 11). The vertical dotted lines give the median value 
for each stage. The median dust luminosity is lowest for the 'weakly interacting' Stage 2 galaxies and increases by more than an order of 
magnitude for the 'strongly interacting' Stage 4 galaxies. The 15-25 K dust temperature is noticeably higher in the Stage 4 galaxies. 
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Fig. 20. — Histograms (top) and cumulative distributions (bottom) of SFR (left) and sSFR (right) derived by MAGPHYS for each 
interaction stage and the control sample shown with the same color and linestyles as Figure 19. SFR shows an increase between non- 
interacting through moderately interacting galaxies and the Stage 4, 'strongly interacting' systems, an evolution similar to that of dust 
luminosity. In contrast, the sSFR distributions are very similar, as is supported by the results of the K-S tests. 
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Fig. 21. — MAGPHYS fits witli tlie 6 different data sets for MlOl, as a representative example of tlie set of fits done for eacli galaxy, 
showing the increasingly constrained parameters (from lower right to top left) as more data are used. UV photometry (present in the left 
panels) is crucial in constraining sSFR, while SPIRE data (present in the top two panels) is essential for the determination of cold dust 
temperature and dust mass. The lower right panel demonstrates the limited constraints that ground based photometry alone can provide. 




Fig. 22. — Parameter confidence intervals increase with the omis- 
sion of data. Points indicate the median size of the 68% confidence 
interval in the probability distribution function (PDF) for each pa- 
rameter for the whole set of galaxies relative to the range when all 
of the data are used in the fit. Using all available data, the median 
68% ranges are: 6.3x10^ L© (dust luminosity), 1.8x10^ M© (dust 
mass), 0.80 K (cold dust temperature), 4.4 K (warm dust tem- 
perature), 2.1 xlO^ Mq (stellar mass), 0.047 M© yr-^ (SFR), and 
1.3 xlO"-*^-*^ yr""*^ (sSFR). UV data are crucial for the determination 
of sSFR, whereas dust luminosity, mass, and 15-25 K temperature 
are best constrained by SPIRE data with further constraints ap- 
plied by photometry from 30-200 /im. 

of these parameters. Figure 21 shows a representative 
example: the best-fit SEDs for all six fits for MlOl as 
well as the accompanying PDFs for the parameters of 
interest. 

As expected, the SFR is significantly constrained by 
UV data. This is demonstrated in Figure 22, where the 
median 68% confidence interval size is a factor of ~4 
larger in the absence of UV photometry. While stel- 
lar mass is fairly well constrained by the UBV-IRAC 
data alone, the absence of UV data tends to result in 
younger stellar population templates being selected by 
MAGPHYS. This effect can be seen in the differences 
in the UV slope and the strength of the Lyman and 
Balmer breaks in the various panels of Figure 21. When a 
younger stellar population template is selected, a smaller 
fraction of the stellar emission is assumed to originate 
from late-type stars, resulting in a tendency to estimate 
the stellar mass ~ 10 — 20% lower than when all the 
data are used. This can be seen in Figure 23 where we 
plot the median fractional difference in the value of the 
galaxy properties for fits with incomplete data sets. UV 
photometry constrains both the SFR and stellar mass, 
and it is also the most important wavelength regime to 
constrain the sSFR. The absence of UV data also tends 
to result in an over-estimation of the SFR resulting in an 
estimate of the sSFR ~40% higher than in fits using all 
available data. 

Herschel data are particularly crucial in constraining 
the cold dust temperature. PACS data typically outline 
the peak of the IR emission, but in cases of the coldest 
dust temperatures, PACS 170 fim is typically too inde- 
terminate and it is only in conjunction with the SPIRE 
250 jam data that the cold dust temperature is reason- 
ably constrained. In contrast, as expected, warm dust 
temperature is not well constrained by the SPIRE data. 
Dust luminosity is typically well-estimated with a com- 



FlG. 23. — Median fractional difference for each parameter rela- 
tive to its value that parameter when all of the data are used in 
the fit. The absence of UV data results in sSFR ~40% higher than 
in fits using all the data. Dust mass, made up primarily of cold 
dust, is very sensitive to changes in the cold dust temperature. 

bination of IRAS and PACS data, only becoming about 
a factor of ~2 more uncertain in absence of SPIRE data. 

SPIRE observations are crucial for constraining the 
dust mass, whose 68% confidence interval would be at 
least a factor of 3 larger without SPIRE data. When 
MIR data at wavelengths A > 30 jam are likewise absent, 
the dust mass becomes almost completely unconstrained 
as the SED contains little information about the dust 
emission. The dust mass estimate is the most sensitive 
to the absence of specific datasets. Interestingly, the dust 
mass estimated by MAGPHYS is ~60% higher when 
SPIRE data are excluded as compared with when all 
data are used; however, when photometry at wavelengths 
A > 30 /im is excluded, the dust mass is estimated ~20% 
lower than when all data are used. The over-estimates in 
dust mass are correlated with the under-estimates in the 
cold dust temperature. Since, cold dust tends to make 
up the bulk of the dust mass and because cold dust mass 
varies as (assuming it is modeled as a /3 = 2 mod- 
ified blackbody), a 5% change in the dust temperature 
results in a 30-40% difference in the dust mass. Aniano 
et al. (2012) recently showed similar results for NGC 628 
and NGC 6946, where fits undertaken only with data at 
wavelengths A < 170 /im tend to overpredict the emission 
in the SPIRE bands and the associated cold dust mass. 
When all data at wavelengths A > 30 /im are omitted, 
MAGPHYS relies primarily on more common UV-NIR 
dominated sources. 

5.4. Comparison of the SFR Derived from the SED to 
Monochromatic and Broad-hand Relations 

Because the SFR is such a key physical parameter, 
many statistical heuristic attempts have been developed 
to estimate it from simple observables such as the UV 
flux (as measured by GALEX) or Ha, the 24 /im flux, 
the total IR flux, and combinations of UV and MIR fluxes 
(e.g., Calzetti 2012). We compared the SFR derived by 
MAGPHYS to the expected SFR based on a wide range 
of relations (Figure 24): from the FUV relations of Salim 
et al (2007) and Treyer et al (2010), from the MIPS 
24 /im emission relations of Rieke et al (2009) and Zhu 
et al (2008, as given in Calzetti et al 2010), from the 
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Fig. 24. — Comparison of the SED-derived SFR to those from 
various SFR proportionahty relations from the hterature. The sohd 
hne shows where the points should lie if the relations agree with the 
SED-derived values, which is the SFR averaged over 100 Myr. The 
dotted, dashed, and dot-dashed lines shows SFR averaged over 20, 
30, and 50 Myr, respectively. SFR is modeled as an exponential 
burst with a peak value such that the average over 100 Myr is given 
by the MAGPHYS SFR and decay times of 20, 30, and 50 Myr. 
While the relations depending solely or partially on IR emission 
agree well with the MAGPHYS SFR, the FUV relations are typi- 
cally high and agree better with a SFR modeled as an exponential 
burst with a width of 30-50 Myr. 

relations combining UV and 24 /am emission of Leroy et 
al. (2008) and Kennicutt et al (2009), and from the total 
IR luminosity relation of Kennicutt (1998). We correct 
the FUV emission from dust using the prescription given 
in Salim et al. (2007) and use the dust luminosity derived 
by MAGPHYS for the total IR luminosity. We assume a 
Chabrier (2003) initial mass function (IMF) and there- 
fore have applied correction factors of 1.06 for relations 
that assume a Kroupa (2001) IMF and 1.60 for relations 
that assume a Salpeter (1955) IMF, following Calzetti 
(2012) and Schiminovich et al. 2007. Figure 24 shows 
the resulting values plotted against the SFR determined 
by MAGPHYS. 

One trend is quickly apparent: the SFRs estimated 
from UV photometry alone tend to be high compared 
with the SED-derived SFR, whereas the relations based 
solely or partially on IR photometry agree fairly well, 
at least for SFR greater than ~ 0.1 Mq yr~^. Median 
differences are of 0.7-0.9 M© yr~^ for the FUV relations. 
Interestingly, our most active galaxy, NGC 3690/IC694, 
shows the inverse trend, as does a prototype starburst 
M82. This may indicate that the correction for dust is 
insufficient for these systems. The outlier of NGC 3190 
(Fig. 24) is due to the low SFR associated with the 
best MAGPHYS model for this galaxy that significantly 
underestimated its UV emission, likely due to the edge- 
on geometry of the system. 

The apparent over-estimation of the SFR by the FUV 
relations is rooted in the time over which the SFR is es- 
timated in MAGPHYS. FUV emission is dominated by 
star formation within the past 50 Myr, although mid- 
to-late B stars can also contribute a significant fraction 
(Calzetti 2012). In contrast, the IR relations typically 
assume that a fraction of the stellar light is absorbed 
to heat dust, and as a result, while the youngest and 
hottest stars dominate the heating of hot dust, the accu- 
mulation of low-mass stars contribute significantly to the 



heating of the more diffuse dust. Hence, the IR relations 
represent star formation over a longer timescale. The 
SFRs derived by MAGPHYS are averaged over the last 
100 Myr, which is more consistent with the timescales 
associated with the IR-dependent SFR relations. Figure 
24 also shows the SFRs averaged over shorter time pe- 
riods with SFR modeled as an exponential decay with 
a peak value such that the average over 100 Myr is the 
MAGPHYS SFR value and decay times of 20, 30, and 50 
Myr, respectively. The SFRs estimated from the FUV 
emission agree better with an exponential decay star for- 
mation history with a width between 30 and 50 Myr, 
which is the expected timescale of a starburst episode. 

6. SUMMARY 

We modeled the FUV-FIR SEDs of fourteen groups 
of thirty-one interacting galaxies, typically with 15—25 
flux points, to determine the most probable evolution 
of dust luminosity, star formation rate, specific star for- 
mation rate, dust mass, stellar mass, and dust tempera- 
ture. The systems were classified as either weakly, mod- 
erately, or strongly interacting (Stages 2-4 in the Dopita 
et al. (2002) scheme). The broad similarities in SED 
shape between different stages emphasize one key conclu- 
sion from this study: as the interaction progresses, and 
even as bursts of star formation may occur, the changes 
are most clearly seen not in the distribution of energy 
broadly but in minor and subtle changes to the SED 
shapes. Bulk SED properties change little, and only 
gradually, in typical interactions. Strongly interacting 
galaxies typically have SEDs characterized by stronger 
MIR emission relative to both their NIR and FIR emis- 
sion and more UV emission relative to their NIR emis- 
sion. 

There are marginally statistically differences (as deter- 
mined by a K-S Test) in the derived galaxy properties: 
dust luminosity and mass, SFR, and cold dust temper- 
ature increase from Stage 3 to Stage 4, SFR increase 
from Stage 2 to Stage 4, and dust luminosity, SFR, and 
cold dust temperature increases from the non-interacting 
galaxies to the Stage 4 galaxies. In contrast, the sSFR 
does not show variations with interaction stage. The rel- 
ative constancy of the sSFR between the different stages 
suggests that this lack of evolution is not due to uncer- 
tainty in stage classification or in the association of inter- 
action stage and progress along the interaction. Rather, 
our set of interacting galaxies shows no clear evidence for 
a burst of star formation prompted by the interactions 
or that such effects occur on timescales such that we see 
enhancements in both stellar mass and SFR, leaving the 
sSFR relatively stable. This suggests a need to be cir- 
cumspect about this canonical activity during the early 
stages of galaxy mergers. 

Different wavelengths have different effects in con- 
straining galaxy parameters in the MAGPHYS SED 
analysis. UV data are important in accurately deter- 
mining stellar population age and that they are impor- 
tant contributors to determining accurately the (specific) 
star formation rate. The stellar mass is primarily deter- 
mined by UBV-IRAC data. SPIRE data are crucial in 
determining the dust mass; in its absence the cold dust 
temperature tends to be underestimated, because the lo- 
cation of the peak is much less constrained without data 
longward of 170 fim. Cold dust mass, which tends to 
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make up the bulk of the dust mass, goes as (assum- 
ing /3=2), so the change in the dust mass is large for even 
a small underestimate in the dust temperature. The pos- 
sible contributions to the SED from AGN are modest for 
this sample and do not affect our conclusions. 

The SFRs derived by MAGPHYS agree reasonably 
well with simple relations based solely or partially on 
IR photometry. Relations based on corrected FUV emis- 
sion tended to overestimate SFR compared to the SED- 
derived SFR, which is averaged over 100 Myr. The SFR 
estimated from FUV can best be understood if it repre- 
sents an exponential decay star formation history with a 
width of 30-50 Myr. 

The complete SIGS sample will bring a significant in- 
crease in statistical power in determining galaxy property 
trends. In addition, testing the accuracy of MAGPHYS 
against simulations of interacting galaxies will help im- 
prove the diagnostic power of SEDs. In a future paper, 
we will examine what kind of simulated interactions and 
their parameters best reproduce observed systems and 
their SEDs. Further, we will test how well MAGPHYS 
recovers galaxy parameters as a function of their inter- 
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APPENDIX 
CLASSIFICATION SCHEME 

The classification methodology used to determine the interaction stage of each system is based on the Dopita et al. 
(2002) classification scheme. Stage 1 galaxies are non-interacting and Stage 5 galaxies are post-merger or coalescence 
systems. Stages 2-4 are weakly, moderately, and strongly interacting systems, based on their degree of morphological 
distortion. In Figure 25, we show a representative system from our sample. The Stage 2 galaxies show little distortion 
and the galaxies are typically well separated. The Stage 3 galaxies have a range of proximity and show some degree of 
distortion and the Stage 4 galaxies show significant morphological distortion and are typically close even in projection. 
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Fig. 25. — Representative examples of Stage 2-4 (left-right) galaxies, showing increasing morphological distortion. 



PROCESSING CONCERNS 
PACS and High Pass Filters 

In the course of processing the PACS photometry, we initially used the PhotProject pipeline scripts to make mo- 
saics. When we compared the measured fiuxes to values obtained by the MIPS and IRAS instruments at the same 
wavelengths, we found some were more than a factor of 2 too low, with MlOl having the greatest disagreement. This 
disagreement is the result of the range of angular extents of the galaxies. The PhotProject pipeline includes a high 
pass filtering algorithm that removes a significant portion of the extended emission in the field and therefore affects 
nearby galaxies with more extended emission much more egregiously. As a result, we found it essential for nearby 
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galaxies to use the Madmap pipeline. 

Swift UVOT and Coincidence Losses 

As we noted in §3.3.2, the UVOT telescope is vulnerable to coincidence losses. When multiple photons arrive in the 
same pixel within the same frame, only a single photon is counted. Losses become greater than 1% when the count 
rate exceeds 0.007 counts per second per pixel. While corrections for these losses and the uncertainties involved is 
well determined for point sources (Poole et al. 2008, Kuin & Rosen 2008), there is a lack of similar understanding for 
extended sources. Although some of the sample galaxies, notably the starburst galaxies M82 and Arp 299, have count 
rates high enough over most of their surface that would result in significant coincidence losses, even for the rest of 
the sample have regions of high fluxes, typically limited to the nuclear area and star forming regions. We therefore 
opted to analyze only those galaxies lacking GALEX photometry as well as an additional test case, as described in 
§3.3.2. Fortunately for the two galaxies having UVOT, but lacking GALEX, only the nucleus and a few small regions 
are bright enough to require corrections. We excluded the regions with high count rates and then measured them 
independently as point sources, correcting them for coincidence losses using the method of Poole et al. (2008). These 
corrections always account for less than 1.5% of the count rate. Based on our test case of NGC 3424, the agreement 
between the GALEX and UVOT photometry is excellent. 

NOTES ON INDIVIDUAL GALAXIES 

Group 1: 

NGC 2976: A small galaxy on the outskirts of the M81 group, NGC 2976 has two active, infrared bright regions at 
the edges of its disk, whose presence is almost completely missed in the UV and NIR. 

NGC 3031: Also known as M81, this galaxy is in our closest system and has very clear spiral arms and filamentary 
structure seen in all the images across our whole filter range. Its nucleus is a LINER. The best-fit decomposition of the 
IR spectrum performed with DECOMPIR has a 4% AGN contribution to the total IR emission and a 16% contribution 
in the 8—35 fim range. 

NGC 3034: M82 is one of the nearest starburst galaxies and has a strong galactic outfiow perpendicular to its stellar 
disk, which is visible in both the UV and Herschel images. Its IRAC 5.8 /im, IRAC 8.0 /im, and MIPS 24 /im images 
are saturated. Its SED is among the worst fit by MAGPHYS with the current set of models. This may be due in part 
to the high obscuration - its fit has a<Ty>=2.9-or perhaps to the inclusion of some emission from the outflow, 
which MAGPHYS does not model. We measured the contribution of the outflow above and below the disk of the 
galaxy within the aperture and found that it contributes ~20-30% of the UV emission and ~ 10-20% of the emission at 
wavelengths greater than 150 /im, which does not fully account for the discrepancy in the UV bands. We selected the 
3.6 /im-derived aperture to minimize the contribution of the emission from the outflow. The outflow's UV emission is 
primarily scattered light from the disk, while its IR emission is from the dust in the outflow. Since MAGPHYS cannot 
model such a feature, its inclusion would tend to bias the best-flt model. 

NGC 3077: Another small galaxy in the M81 group, NGC 3077 lies behind a nearby bright star that prevented it 
from being observed in the UV by GALEX and UVOT. Its MIR-FIR images shows evidence of tidal stripping in the 
asymmetric structure and southwest lobe. 
Group 2: 

NGC 3185: This Seyfert galaxy shows a circum-galactic ring of star-forming material in both UV and IR images. 
The DECOMPIR decomposition requires a 3% and a 12% contribution of the AGN to the total IR luminosity and 
8—35 /im luminosity respectively. Its disk shows up in the FIR image as two bright regions at the opposite ends of 
the galaxy. It is the most distant member of the triplet it forms with NGC 3190 and NGC 3187 on the sky, but it 
and NGC 3190 have much closer recessional velocities than NGC 3187 (1217 km s"^ and 1271 km s ^ vs 1581 km 
s~^). Group 2 is the only compact group in our sample, although the SIGS sample has several others. Tzanavaris 
et al. (2010) found a bi- modality for the sSFR distributions of compact groups depending on the slope of the IRAC 
photometry between 4.5 /imand 8.0 /im. NGC 3185 has a negative IRAC slope and its sSFR places it in the relatively 
quiescent population as expected by Tzanavaris et al. (2010). 

NGC 3187: This galaxy has a pair of tidally elongated arms, which are best detected in the UV. This galaxy did not 
have archival IRAS fluxes. NGC 3187 has a positive IRAC slope and its sSFR places on the edge of the distribution 
of galaxies with positive slopes, argued by Tzanavaris et al. (2010) to indicate active star formation. 

NGC 3190: The most massive of the three galaxies in this group, NGC 3190 is a nearly edge on LINER and has a 
dust lane, that appears prominently in absorption in UV and is correspondingly bright in IR emission. Its SED shows 
particularly low UV relative to its visible emission, but presumably the geometry of this dust lane explains the poor 
flt to the UV and the disagreement in the estimates of its SFR. The DECOMPIR best-flt requires no contribution 
from an AGN. NGC 3190 has a negative IRAC slope, but the low SFR value derived by MAGPHYS places it in the 
quiescent category of Tzanavaris et al. (2010). 
Group 3: 

NGC 3226: As an elliptical galaxy, NGC 3226's emission is dominated by its stars, although there appears to be 
a faint tidal feature in the 8 /im emission directed roughly to the north. NGC 3226 is particularly faint in the MIR 
relative to its NIR emission and only has upper limits in the MIR-FIR from IRAS. It has a LINER nucleus, but the 
decomposition of its IR spectrum indicates no signiflcant contribution form an AGN. 

NGC 3227: Along with its smaller companion NGC 3226, NGC 3227 was not observed by GALEX, so we use UVOT 
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data instead. Its IRS spectrum shows [Ne v] emission consistent with a Seyfert nucleus. Its DECOMPIR fits are the 
worst of the nine galaxies, but have the highest AGN contribution of 15—25% to the total IR luminosity and 45—60% 
of the MIR luminosity. 
Group 4: 

NGC 3395/3396: The apertures of this pair of galaxies overlapped sufficiently to make determination of the emission 
belonging to each galaxy problematic. We therefore opted to treat the system as a combined system. This pair of 
galaxies is distantly associated with Group 5. This group does not have UBV photometry. The SED shows little 
attenuation in the UV and strong 60 /im emission relative to the 100 fim emission. 
Group 5: 

NGC 3424: This edge-on galaxy's central region becomes increasingly bright relative to its disk with increasing 
wavelength. Similarly to NGC 3190, the other nearly-edge on galaxy, its NUV-FUV slope is quite steep and not 
particularly well fit by MAGPHYS. We used this galaxy as a test case for analyzing UVOT photometry and found 
good agreement with the GALEX photometry. It lacks UBV photometry. 

NGC 3430: Seen nearly face-on, this galaxy has a large and extended UV disk. This system provides a nice example 
of a system fairly early in the interaction sequence and hence fairly undisturbed morphologically. The IR peak of 
its SED is not very well constrained as this galaxy was not observed by PACS and did not have MIPS 70 fim and 
160 /im fluxes available. 
Group 6: 

NGC 3448: This is the larger member of this pair of dwarf galaxies. A bridge of emission extends from NGC 3448 
in the direction of UGC 6016, seen most prominently in the NUV. 

UGC 6016: This dwarf galaxy is very faint in the IR, but has a large but diffuse UV envelope. Due to its lack of 
significant detections in several of our filters and in ancillary IRAS^ we cannot constrain the galaxy parameters very 
tightly. Its distance is also the most uncertain of the sample as it did not have a distance modulus or a recessional 
velocity in the PSCz catalog. 
Group 7: 

NGC 3690/IC 694: Also known as Arp 299, this system is the most active in our sample, with the highest amount 
of star formation and is our only LIRG, showing a corresponding large amount of UV attenuation. Its 8 /im image is 
saturated, and it does not have ancillary UBV fluxes. 
Group 8: 

NGC 3786/88: The galaxies in this pair have very similar UV— NIR fluxes, but NGC 3788, which is more edge-on, 
has higher fluxes in the Herschel bands. However, their best-fit SEDs have similar infrared luminosities, likely due to 
the relatively high contribution of the warm dust in the model for NGC 3786 compared to NGC 3788. Neither galaxy 
has ancillary IRAS or UBV fluxes. NGC 3786 is a Seyfert galaxy, with a bright nuclear region showing [Ne v] emission 
in its IRS spectrum and a partial ring of star-forming regions. Its decomposition requires an AGN contribution of 7% 
and 26% to the total IR luminosity and MIR luminosity respectively. 
Group 9: 

NGC 4038/4039: Also known as the Antennae, this pair of galaxies is one of our most evolved systems and, because 
we cannot separate them, we treat it as a single entity. Its clumpy distribution of star forming regions are clearly 
apparent at 8 /im and in the UV bands; the two nuclei are most clearly seen in 2MASS and IRAC, a reflection of the 
relative PAH strengths to warm dust. 
Group 10: 

NGC 4618: Paired with NGC 4625 in a roughly equal-mass dwarf galaxy interaction, NGC 4618 has an off-center 
nucleus with a single arm curving to the south-west, features seen in all the images. 

NGC 4625: While relatively compact in the infrared, NGC 4625 has a faint, diffuse set of flocculent set of spiral 
arms observed best in the NUV band. 
Group 11: 

NGC 4647: This spiral galaxy is located at the edge of the IRAC field, which complicated the measurement of 
its flux because its aperture, determined on the NUV image, extends past the edge of the IRAC image. We had to 
manually correct for the edge pixels without flux. We also do not have ancillary UBV fluxes. 

NGC 4649: A large elliptical galaxy also known as M60, NGC 4649 is very faint in the mid-infrared and absent in 
the far-infrared, a dramatic contract to its companion. It has the lowest sSFR of the sample. 
Group 12: 

M51A: The Whirlpool Galaxy is the larger galaxy in this well-studied system. It has quite consistent morphology 
across the wavelengths, but with more inter-arm filamentary emission and greater extent in the UV than in the IR. 
The decomposition of its IR spectrum is best fit without an AGN contribution, despite its LINER nucleus. This is 
one of the systems where the flux from the smaller galaxy was subtracted from the aperture of the large galaxy. 

M51B: The smaller companion to the Whirlpool Galaxy, M51B is dominated by early- type stars and has very little 
UV emission. Its FIR emission is confined to its nucleus. Its MIPS 160 /im measurement was a factor of three lower 
than the PACS measurement, and we opted to omit it from the fits. 
Group 13: 

NGC 5394/95: The smaller galaxy of our most distant pair, NGC 5394 shows a beautiful pair of tidal tails, especially 
in the UV. NGC 5395 has a LINER nucleus, which requires only a small AGN contributions of 3% and 12% to the 
galaxy's total IR and MIR luminosities respectively in the DECOMPIR fits. This is the other system where the small 
galaxy flux needed to be subtracted from the large galaxy aperture. 



36 



Lanz et al. 



Group 14: 

MlOl: Another well-studied galaxy, its GALEX^ Spitzer^ and Herschel images show clumpy star forming region 
structures along the spiral arms. Its large size made determining a single unbroken aperture with SExtractor compli- 
cated. As a result, some of the outermost UV emitting regions were not included. 

NGC 5474: A small companion to MlOl, its core is offset to the north from its center. 
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